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^ : Abstract 

, Thermodynamics in the vicinity of a critical endpoint with nonclassical expo- 



a 



nents a, /?, 7, 5, • • • is analyzed in terms of density variables (mole fractions, 
magnetizations, etc.). The shapes of the isothermal binodals or two-phase 
coexistence curves are found at and near the endpoint for symmetric and 



O ■ nonsymmetric situations. The spectator- (or noncriticalj-phase binodal at 

O 

T = Tg is characterized by an exponent {6+l)/6 (2:^ 1.21) with leading correc- 
^ : tions of relative order 1/6 (2^ 0.21), 6>4//35 (^ 0.34) and 1 - {(36)-^ (2^ 0.36); 

in contrast to classical (van der Waals, mean field, • • •) theory, the critical 
endpoint binodal is singular with leading exponent {1 — a)/f3 (^ 2.73) and 
corrections which are elucidated; the remaining, A-line binodals also display 
the 'renormalized exponent,' {1 — a)/(3 but with more singular corrections. 
(The numerical values quoted here pertain to {d = 3)-dimensional-fluid or 
Ising-type systems.) 



^) Current address: Center for Polymer Studies, Center for Computational Science, and Depart- 
ment of Physics, Boston University, Boston, MA 02215 USA 



I. INTRODUCTION AND OVERVIEW 



At a critical point in a fluid (or other Ising-type or n = 1) system two distinct phases, 
say, (3 and 7 become identicah below T = Tc these two phases may coexist for appropriate 
values of the conjugate ordering fleld (or chemical potential, etc.)0 h] above Tc they merge 
into a single phase, say, f3j. If there is some other fleld variable,^ say which may be 
varied without destroying coexistence, the critical point is drawn out into a lambda line, 
T = Tc{g). A typical situation, which lacks any special symmetry, is shown in Fig. |1|. The 
lambda line. A, delimits the phase boundary surface h = hp{g^ T) labeled p, on which (3 and 
7 may coexist. 

Now in many instances when g is varied, say decreased, another quite distinct phase, a, 
will be encountered. In this case the lambda line terminates at a critical endpoint^ which is 
labeled E in Fig. |T|. At E the phases f3 and 7 may undergo criticality in the presence of the 
coexisting noncritical phase a which may be appropriately termed the spectator-phaseE^ 
The surface bounding the spectator-phase in the (g^T^h) or fleld space is labeled a; on it 
a may coexist with phases /?7, /?, or 7; on the triple line, r, where the surface p meets the 
surface a, all three phases a, f3 and 7 may coexist. 

In a previous studyi (to be denoted I), we discussed the shape of the spectator-phase 
boundary surface, g = ga{T^ /i), in the vicinity of the endpoint at T = Tg and, by choice of 
origin, h = = 0. It was found that the surface is singular at E with functions such as 
griT)^ specifying the triple line, and gaiT^; h) displaying nonanalytic behavior described by 
a variety of critical exponents.li@ When, as is normally so, the lambda line is characterized by 
nonclassical exponents a, for the speciflc heat, /?, for the order parameter, 5, for the critical 
isotherm, etc., the spectator-phase boundary exponents can all be expressed! in terms of 
a, /?, and S. Beyond that it was shown that various dimensionless ratios constructed from 
the amplitudes of the phase-boundary singularities should be universal with values also 
determined by the nature of the bulk criticality on the lambda line.ii 

These conclusions were based on a phenomenological description of the thermodynamic 
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potentials (or Gibbs' free energies) G^{g^ T, h) and G^^{g^ T, /i), for the spectator-phase and 
for the coexisting and critical phases, respectively. The former was assumed to have a power 
series expansion in the vicinity of E] the latter embodied a full scaling representation of the 
critical line and its neighborhood.BS 

This formulation neglects the essential singularities expected on the a and p boundaries ;i 
these can, however, be discussedS but play only a negligible quantitative role. Our general 
phenomenological treatment has been checked by an extensive study of a family of spher- 
ical models which exhibit lambda lines and critical endpoints with a range of nonclassical 
exponents (although /? = | in all cases) 

Many experimental examples of critical endpoints are found in multicomponent fluid 
systems. In the simplest example, which we will particularly bear in mind, two chemical 
species, B and C, mix as fluids in all proportions at high temperatures forming the phase /?7. 
At lower temperatures, however, they undergo liquid-liquid phase separation, or demixing, 
producing phases (3 and 7 rich in B and C, respectively. Up to a constant shift, the fleld 
h may then be taken as the chemical potential difference Pb — l^c- As the pressure, p, or 
the total chemical potential, Pb + I^C) either of which we may identify with the fleld is 
reduced, a dilute vapor phase, a, appears. Fig. |1] then represents a characteristic overall 
phase diagram. Now in a typical experiment the temperature T is controlled and may be 
held flxed: corresponding to Fig. |l], the appropriate isothermal phase diagrams in the {g^ h) 
plane then have the character shown in Fig. |^ for (a) T < Tg, (b) T = Tg, and (c) T > T^. 

However, the chemical potentials jiB and /ic, or the flelds h and g^ are normally not 
under direct experimental control or observation; rather, the conjugate densities^ Pb and pc 
(or concentrations of B and C) or, equivalently, the densities 



Pi = -^G{9,T,h) 



d 

and p2 = - -T^Gig, T, h) 
9,T 9g 



(1-1) 

T,h 



are the prime experimental variables. [Note that in the example envisaged with g = Pb + l^c 
one simply has pi = \{pb — Pc) and p2 = \{pb + Pc)-] In the density plane (pi,p2) the 
phase boundaries p and a are represented by two-phase regions bounded by smooth curves. 



the so-called binodals or coexistence curves: see Fig. The aim of this article is to analyze 
in detail and generality the shapes of these isothermal binodal curves in the vicinity of a 
critical endpoint. Specifically, we will elucidate the nature of the leading and subdominant 
singularities that appear in the various binodals labeled + , etc., in Fig. 0. 

It appears from Fig. ^, and detailed analysis bears it out, that the binodal curves for 
T > Tg meet with a common tangent at the endpoints and and at the extended 
triple points fa and ff3j (defined by the intersection of p, the extended phase boundary p, 
with the surface a in the (g^T^h) space: see I). Of principal concern, then, is the way in 
which the binodals depart from linearity. Above Tg one expects analytic binodals but the 
behavior of the curvatures at fa and f/?7 as T ^ Te+ is then of interest. On the other 
hand, atT = one expects singular behavior at E^ and E^. Indeed, Borzi,i stimulated by 
Widom,i discussed the noncnizca/ binodals at the endpoint, namely, B^^ in Fig. |^(b), using 
the simplest possible phenomenological postulate and geometrical arguments (equivalent to 
van der Waals and other classical theories). He concluded that the degree of tangency was 
controlled by a 4/3 power law (in place of a power 2 for a normal analytic tangency). 

Later Klinger,i using a more general phenomenological classical theory, discussed the 
critical endpoint binodals, i3f and analytically: see Fig. ^(b). However, he found no 
evidence of singular behavior. Beyond that, Klinger confirmed the leading 4/3 power in the 
noncritical or spectator binodal and found that the first correction term carries a 5/3 power. 

On general grounds, however, it seems certain that the powers 4/3 and 5/3 must result 
from the reliance on classical theory which entails the critical exponent values a = 0, /? = 
|, and = 3 in place of the appropriate nonclassical values a O.lOe, /? — 0.328, and 
5 = (2 — a)//3— 1 4.78 which characterize the specific heat, coexistence curve, and 
critical isotherm of real bulk, {d = 3)-dimensional fiuids (or other systems in the Ising 
universality class). Indeed, Widom has conjecturedi that in general the 4/3 power should 
become {5 + l)/S. This reduces to Borzi's result when S = 3 but yields an exponent value 
1.207 for real fiuid systems. 

Here we confirm Widom's surmise using the full scaling approach developed in I. Fur- 
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thermore we show that Khnger's correction exponent of 5/3 is replaced, more generaUy, by 
three exponents, namely, (2 — a + /?)/ {2 — a + 64)/ (35 and (3 — 2a — /?)/ (35. Here 6^ is the 
leading correction-to-scaling exponent which has the value ^4 0.54 for [d = 3)-dimensional 
Ising-type systems;^ thus these three exponents have values of about 1.42, 1.55 and 1.57, 
respectively, for bulk fluids. In the classical limit attained via d —> 4:— one has ^4^0 and 
the second exponent reduceJill to 4/3 while the flrst and third yield Klinger's value 5/3. 
However, we also identify further singular exponents that must appear in the expansion of 
the noncritical binodal at the endpoint. 

It transpires, in addition, that, contrary to Klinger's flndings,! the critical binodal is, in 
general, also singular with a leading power [1 — a)/ (3 2.7^ so that the binodal is much 
flatter at the endpoint than classical theory would predict. Here, and below where 
appropriate, we suppose a > as applies to real fluids. The exponent {1 — a)/(3 is, in fact, 
the same as that long known to characterize isothermal binodals passing through a lambda 
point (away from any endpoint): see B^^ and in Fig. ^(a). This behavior which is, of 
course, reconflrmed by our analysis reflects, in turn, the phenomenon of critical exponent 
renormalization.lli The correction terms in the critical endpoint binodal are found to carry 
exponents {l — a + 9k)/(3 with A; = 4, 5, • • •. When one substitutes the classical values a = 
and 9k = — 4) these leading and correction exponents become 2, 3, 4, • • • which are 
consistent with Klinger's results and indicative of a fully analytic critical binodal. 

The results sketched here, and others for the remaining binodals shown in Fig. |, are 
presented in detail in Sec. fT|. However, it is necessary to point out that Figs. ||-3 are 



special in two respects. First, as mentioned, no symmetry with respect to the ordering 
surface p has been supposed: this is quite appropriate for most fluid systems. However, 
as observed in I,3S there are many other physical systems in which the thermodynamic 
potentials are unchanged under reflection in the plane p: one may then take h = on p and 
the symmetry becomes invariance under h <^=^ —h. The conceptually simplest example is an 
elemental ferromagnet, like nickel or iron, where h = H is the magnetic fleld and g = the 
pressure. Other examples are ferroelectrics, antiferromagnets, order-disorder binary alloys. 



and liquid helium through its transition to superfluidityp however, the binodal curves are not 
readily accessible experimentally in some of these cases. The corresponding (g^T^h) phase 
space, the isothermal sections, and the binodal curves for such symmetric critical endpoints 
are illustrated in Figs. ^, ^, and |^. In fact, symmetric critical endpoints are simpler in a 
number of respects and will be analyzed first below. Fundamentally we find that the leading 
singular behavior of the binodals is identical in the symmetric and nonsymmetric cases but 
the correction terms diflFer in character: see Sec. fT|. 

A second special feature embodied in Figs. |l|-3 is the slope of the A line which we char- 
acterize as negative in the sense that if, without loss of generality, we (i) take 

g = h = 0, T = Te, at (1.2) 

and (ii) suppose that the negative g axis lies in the a or spectator-phase (see Figs. 1 and 4) 
then we have! 

Conversely, as illustrated in Fig. |^, one must also consider the case of a positively sloping A 
line with 

As seen in Figs. ||and||, this produces distinct isothermal phase diagrams and new arrange- 
ments of binodal curves: note the additional notation introduced in Fig. |^. 

One might, of course, also wish to consider the borderline cases A^ = 0, oc; we will not 
pursue these but, on the basis of our postulates for the thermodynamic potentials as set out 
below in Sec. the necessary analysis presents no further problems of principle. 

In summary therefore, we will analyze the binodals for four types of critical endpoints 
which, using N for nonsymmetric and S for symmetric, may be labeled NA (Figs. |l]-3), 
NB, SA, and SB (Figs. §6). 

In outline the remainder of the paper is as follows: Our basic scaling postulates for 
the thermodynamic potential G^^{g^T^h) are set out in Sec. ||. They are essentially the 
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same as introduced and discussed critically in I but they have been extended significantly 
as regards the symmetries of the corrections to scaling and of the nonlinear scaling fields; 
the notation also differs in a few details.! The reader prepared to take the postulates on 
tru J may proceed directly to Sec. 01 where the shapes of the birrodals in the various cases 
are discussed in detail without reference to Sec. y. The analytic derivation of the results, 



which is straightforward in principle but a little delicate in practice, is presented in Sec. fV . 
Explicit formulae for the many amplitudes entering the expressions for the various binodals 
in Sec. |T| are also given in Sec. Section summarizes our conclusions briefiy. 



II. THERMODYNAMIC POTENTIALS FOR ENDPOINTS 

This section sets out a complete specification of the thermodynamic potential G{g^T^ h) 
in field variables as needed for the general description of critical endpoints. It is the basis 
for the results described in Sec. |ITI| but need not be read to understand those results. For 
convenience we adopt the critical endpoint as origin for the fields g and h as specified in 
(|L2D, and also put 

t = {T-T,)/T,. (2.1) 

Thus t, and h measure field deviations from the endpoint E at (g^t^h) = (0,0,0). (In 
I the variables g and t were denoted Ag and i.) For any property P{g^T^ h) admitting a 
power series expansion about E (of indefinitely high order but not necessarily convergent) 
we utilize, for brevity, the semisystematic subscript notation 

P{g, T, h) = Pe + Pig + P2t + Psh + P^g^ + 2P,gh + 2Pegt 

+2Prht + Pst^ + Pgh^ + Osig, t, h), (2.2) 

where, here and below, Om{x^y^z) denotes a formal expansion in powers x^y^z^ with j + 
k + I > m. If P is symmetric under h —h one has 

Ps = P5 = P7 = (to order 3). (2.3) 
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Functions satisfying (|2.2| ) and (|Z3|) will be said to be noncritical (as against critical). 

Following I we assume that the thermodynamic potential G^{g^T^ h) for the spectator- 
phase, a, is noncritical. Thus one has, e.g., = ^[d'^G^{g^T^ h) /dhdt]^ and, by virtue of 
(|1 . 1|) , the endpoint densities in the spectator-phase are simply 

pr = -G^ and pr = -Gl (2.4) 

To describe the critical phases, /?, 7 and f3j^ we first introduce, again following I, the 
two relevant nonlinear 'thermal' and 'ordering' scaling fields, i{g^ T, h) and h{g^ T, /i), which 
both vanish on the A line while h also vanishes on the phase boundary p. For the nonlinear 
scaling fields we accept the noncritical expansionslll 

i = t + qoh + qig + q2g^ + q^gt + q^^ + q^gh 

+ qGh^ + q7th + 0s{g,t,h), (2.5) 
h = h + r_it + rog + righ + r2th + r3/i^ + r^g^ 

+ r,gt + ret^ + Os{g,t,h), (2.6) 

which slightly extend those in 1(4.7), (4,8). It should also be mentioned at this point that 
pressure-mixing terms, which have been discovered recently in connection with the Yang- 
Yang anomaly in fiuid systems,0'lll are not considered here.0 
In the symmetric case one has, to order 3, 

Qo = Q5 = Qr = and rj = for j = —1, 0, 3-6. (2.7) 

Asymptotically, the A line may thus be described by 

gx{T) = Agt + A,2t' + O(t^), hx{T) = A^t + Ah2t^ + O(t^), (2.8) 

where one finds 

. .l^JHr^, A, = ri^ilTzl, (2.9) 

<7i - 90^0 qi - qoTo 



8 



with similar expressions for A^2, etc. In accord with ( |1.3|) and ( pTiP , we assume does not 
vanish or diverge. In the symmetric case one has = —1/qi and = Ah2 = • • • = 0, so 
that Ao = gi - qoTo 7^ 0. 

Then we need the one relevant scaled variable 

y{g, t, h) = Uh/\i\^, with A = /35 = /3 + 7>l, (2.10) 

where the exponent relations and inequality are standard. In I we took U = U{g^t^ h) as a 
noncritical function; however, with no loss of generality we may take [/ as a positive constant 
since any dependence on and h can be absorbed into h. Beyond y we need the many 
irrelevant scaled variables 

yk{g,t,h) = Uk{g,t,h)\i\'\ 

ek+i>Ok>0, fc = 4, 5,.... (2.11) 

We assume that the associated irrelevant amplitudes Uk are noncriticallll with 

Uk{g, t, -h) = {-fUk{g, t, h) in case S. (2.12) 

Now we can write the thermodynamic potential for the critical phase as 

G^\g, T, h) = G\g, T, h) - g|tT""W^±(?/, 1/4, V,, ■■■), (2.13) 

where the background G^{g^ T, h) and the positive amplitude Q{g^T^ h) are noncritical while 
the subscript ± refers to t ^ 0. Physically, from the relation of a to the specific heat we 
have 2 — a > 1 but we further suppose 

{2-a)/A = {S + l)/S>l, (2.14) 

as is generally valid both classically and nonclassically. For concreteness and simplicity we 
will, in addition, focus on a > (as appropriate for bulk fiuids, etc.). 

We also assume, acknowledging the symmetry of the standard universality classes, that 
the scaling function W±{y^ 2/4, ?/5, • • •) is both universal and invariant under change of sign of 
the odd arguments ^, ^5, 7/7, • • •. Beyond that we have the expansion 
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W^{y, 1/4, ?/5, ■ ■ •) = Wl{y) + yM\y) + y^w!i\y) + • ■ • 

= Y.^^^y)y^'^^^ (2-15) 

in terms of the irrelevant scaled variables 2/4, 2/5, • • where for brevity we have introduced 
the multi-index 

K = 0, (4), (5), • ■ • , (4, 4), (4, 5), ■ ■ • , (4, 4, 4), • • ■ , (2.16) 

and the associated conventions 

l/° = l, - (2.17) 

We also say k = [(i, jr ' ' ^ ^)] is odd or even according as the sum i + j + • • • + n is odd or 
even. Then with an obvious extension of notation, the symmetry of W±{y^ • • •) requires 

W^i-y) = H'^W^iy). (2.18) 
For small y and t > we can then write the further expansions 

W^{y) = ly^o + y'W^2 + y'W^A + • • • , for A. even, 

= + y^W^^ + y^W!^^ + • • • , for odd. (2.19) 

These series may, in general, be normalized via 

W"+2 = ^+0 = 1 ('^ even) or W^-^ = 1 {k odd), (2.20) 

which serve to fix the nonuniversal metrical amplitudes Q, U and Uk^e^ etc. 

Note, however, that in setting W^q = = +1 an appeal to thermodynamic convexity,i 
together with Q > and a > 0, is entailed: see Ref. 17 where the consequences of the 
necessary convexity of the basic thermodynamic potentials is discussed both for the scaling 
functions and, more generally, for critical endpoints, thereby extending Schreinemakers' 

ruies.liiii 
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For i < the existence of the first-order transition leads to \y\ factors in the expansions 
so that one has 

W^iy) = [W'^o + \y\W^i + y'W^ + \yfW^, + • • ■]aK{y), (2.21) 
where the special signum function is defined by 

^Hi{y) = 1 for even, 

= sgn{y) for k odd. (2.22) 

Convexity with U > then shows that W^^ and W^2 i^^st both be positive: see Ref. 17. 

For large arguments, \y\ oc, the individual scaling functions W^{y) and W'^{y) must 
satisfy stringent matching conditions to ensure the analyticity of G^^ {g^T^h) through the 
surface t = for all /i 7^ 0. These often overlooked conditions may be written 



l + Y.w'^{±\y\)- 



-l/A 



1=1 



<jM, (2.23) 



where the multiexponent 9[k] is defined by 

^[0] =0, e[{i,jr",n)]=e, + 9j + --- + 9r,, (2.24) 



with i, j, • • •, n > 4. By virtue of the normalizations (|2.20|) the numerical amplitudes 
W^-, and wj^ should all be universal (as should the exponents, a, /?, 5, ^4, ^5, • • •). 

Beyond that, as shown in Ref. 17, convexity dictates that and W2 must be positive while 
{w^y/w2 must be bounded above. The sign of is not determined by convexity alone but 
must, in general, be negative: see Ref. 17. This plays an important role in determining 
allowable density diagrams. 

Finally, from ( p.l3|) we note that the critical endpoint densities are 



p^,'' = -Gl and p^^ = -G?; (2.25) 

see Fig. 3(b). 
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To close this section we recall from I that the phase boundary a oi g = ga{T^ /i), follows 
by equating the two expressions G = G^{g^T^ h) and G = G^^{g^T^ h). Consequently, it is 
useful to define the thermodynamic potential difference 

D{g, r, h) = G^ig, T, h) - G\g, T, h), (2.26) 



which is noncritical by virtue of the definition of G^ in ( |2.13|) . By our conventions the 



negative g axis, i.e. t = h = 0^ g < 0^ lies in the a phase (see Figs. [T| and ||): this implies 
Di > 0. The phase boundary p and its extension p above Tc{g)^ is given by h{g^ T, h) = 0. 
As in 1(5.4), we will assume that the A-line is not tangent to the triple line r at E. The 
densities pi and p2 on the boundaries a and p then follow from (|Ll|) and, by eliminating 
g and h at fixed T, the various isothermal binodals can be computed as expansions about 
E or about A: see Figs. |^ and ||. We postpone the details until Sec. |I^ and turn next to 
describing the results. 



III. THE ENDPOINT BINODALS AND THEIR INTER-RELATIONS 

In this section we describe the results of our analysis of the possible shapes of the various 
binodal curves and their inter-relations with one another as illustrated in Figs. || and |^. 
After some preliminaries describing the "rectification" of the binodals, we consider first 
the behavior near the A line: this entails only the free energy G^^{g^T^h) and, inasfar as 
the corrections to scaling are involved, extends previous knowledge somewhat. Then the 
binodals at the critical endpoint temperature T = Tg are described: these are, perhaps, 
of most interest. The binodals associated with the a surface above Tg are discussed next. 
These are analytic but their slopes and curvatures display critical singularities as T ^ Te+. 
Finally, the binodals associated with the three-phase triangle below Tg are considered. 
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A. Rectification of the Binodals 



We approach the description of the binodal curves by supposing that at fixed T one 
may observe the densities (pi, P2) of various pairs of coexisting phases. In binary fluid mix- 
tures pi and P2 might correspond directly to the number densities of the two species, B 
and C. In ternary mixtures^ however, observations would normally be conducted at fixed 
temperature and pressure and varying composition. Then pi and p2 would each represent 
convenient linear combinations of the number densities of the three species, say. A, B, and C 
as represented typically in a triangle diagram. (Our analysis also applies to observations of 
quaternary mixtures if sections of the thermodynamic space corresponding to constant tem- 
perature, pressure, and a third field (or combination of chemical potentials) are constructed; 
however, experiments are not normally conducted that way and some further analysis would 
be needed to describe, say, a section at constant T, and ^3.) 

We suppose next that the critical endpoint temperature Tg itself can be determined with 
reasonable precision so that the variable t = {T — Tq)/Tq of ( p.l|) is well defined. Then the 
densities (Pi^,P2^) = and (Pi^,P2^) = of the spectator and critical phases at the 
endpoint can be found: see Figs, ^(b) and |§(b). These define an axis of slope 

= Api/Ap2 = (P^^ - pT)/{P2 - pT)- (3.1) 

A natural second axis is found by noting that according to classical theory§ the critical 
binodals i3f and BJ have a well defined common tangent at E^ of slope {dpi/ dp2Y (3 = 
1/Lp, say. This is confirmed by our more general analyses which, indeed, predict that the 
binodals are fiatter at E^ which eases the practical determination of Lp. (Note that it proves 
convenient to define Lp reciprocally with respect to L(j\ see below .0) 

To describe the various binodals near the endpoint it is then natural to adopt new density 
variables, m and m, which are linearly related to pi and p2 but utilize E^ as origin and are 
oriented along the axes just specified: see Figs. |^(b) and|^(b). Henceforth, therefore, we will 
utilize the rectified density variables 
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m = pi-pY-L^{p2-pT). (3-2) 
m = p2-p^'-Lp(pi-p^'). (3.3) 

Furthermore, without loss of generahtyH we assume that the only pure phase located within 
the quadrant m > 0^ m > 0^ at T = is the (3 phase. Then, as iUustrated in Figs, ^(b) 
and ^(b), the a phase at T = Tg is restricted to m < and only the 7 phase lies in the 
quadrant m < 0, m > 0. 

The notation m and fa is suggested by the magnetic case in which m, the magnetization, 
is the primary order parameter discontinuous across p and coupling to the ordering field 
/i, while m is a secondary or subdominant order parameter conjugate to g. Note, indeed, 
that for symmetric endpoints we have@ L^r = I/p = so that if one shifts the definitions of 
the densities in a natural way to yield an origin = = one simply has m = pi and 
fh = P2: see Fig. 

B. Lambda Line and Associated Binodals 

We note first (that within the postulates of Sec. HD the densities on the A line are 
noncritical functions of T so that we have 

mx{T) = Mit + M2t^ + ---, (3.4) 
mx{T) = Mit + M2t^ + • • • . (3.5) 

For a symmetric endpoint all the Mj vanish identically. Beyond that, the coefficients Mj 
and Mj are not restricted in magnitude or sign although, of course, the A line itself cannot 
extend beyond the endpoint. Thus one must, here, have t < in case A and t > in case 
B. 

Next notice that the binodals B^^ for T < Tg, and — see Fig. 1^/3^^ for T = Te, and 
By"^ for T > can all be treated together since by our postulates all of these binodals 
depend only on the free energy of the critical phase. Furthermore, inasfar as they are 
not truncated by the spectator-phase, they must all share the same singularities and vary 
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uniformly with T. It is also convenient to describe the binodals with the aid of a parameter 
s > (related to which vanishes on the A line and increases into the /? and 7 phases: 
coexisting phases correspond to the same value of s. 

In the symmetric case^ S, the binodals associated with the A line or p boundary may 
then be specified by 

m± = ±Bs [1 + hs'^'^ + hs"^ + 

+ . . . + fesZ+C'^/^) + •■•], (3.6) 
fh = rhx{T) + [1 + 04/^/^ + ais^/^ + • ■ • + + ■ • 

+ks^/^ [1 + hs"'^ + ■ • ■ + hs^'^ + ■•■]. (3.7) 

In ( |3.6| ) the general correction term has the form 6n(^)s^*-"-' where n = [n^] is a multi- index 
with > and the exponents here and in ( p.7| ) have the form 

/3C(n)= no + 5^712,^2,, (3.8) 

J>2 

/?C(n) = no + J][n2,^2, + n2,+i(A + O^j+r)]. (3.9) 

The appearance of the exponent A = /?5 is due to the symmetry which acts to suppress the 
odd irrelevant variables. 

The correction amplitudes a4(t), ai(t), • • •, 64(^)5 ••• are noncritical but, generally, of 
indeterminate sign. However, the noncritical amplitude B{t) = + + • • • is positive 
with our conventions and the signs ± correspond to the (3 and 7 phases, respectively. The 
amplitudes A{t) = Ae + + • • • and K{t) = + + • • • are also noncritical. For a > 0, 
as we may assume here, the amplitude A must be negative in case A while it is positive 
in case B. For a < the amplitude K would have to have matching signs but that is not 
demanded for a > 0. Explicit expressions for A^^ B^^ etc., are given in ( [4.26|) and (|4.27| ). 

It is clear by symmetry that the (m, m) tielines connecting coexisting phase points are 
all "horizontal," that is, parallel to the m axis (m = 0): see Fig. |^. Similarly, the diameter of 
the p binodals, defined as the locus of midpoints of the tielines, is given simply by nidiam = 0, 
rhdiam >rhx{T)>0. 
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The symmetric A binodals may finally be expressed directly in terms of m as a variable 
by solving (|3.6| ) for s and substituting in (|3?7|) . With x = \m/B\ this yields 

m = rhx{T) + [l + a^x^'/^ + aix^/^ + • • •] 

+ Kx^/^ [1 + ^4^^^"/^ + hx^^^ + •••], (3.10) 

where (24 — (24 — (1 — a)h/^/ (3 and so on. The term in A provides the dominant behavior 
(when a> Q) with {1 — a)/ (3 2.73 for Ising (i = 3 as quoted in the Introduction. However, 
the term in K provides strongly competing corrections of relative order \m\^^^: note that 
1/(3 3.05. The higher order correction terms run through all powers of x with exponents 
of the form C(ni) + ({^2) + ••• + Ci^i)- 

The nonsymmetric^ N, binodals associated with A line or p surface — see Fig. |^(a) — 
may be described similarly. In terms of the parameter s we find 

m± = mx{T) ±Bs[l + hj''^ + h^s^'^ + • • • ± h^s^''^ ± • • •] 
+^^(1--)/^ [1 + aj^'^ + ars"^ + • • • ± a,s'^'^ + • • •] 

+Ks^'^ + + + , (3.11) 
m± = fhx{T) + is^^-")/^ [1 + aj''^ + ai5^/^ + • • • ± a^s^^'f^ ± • • •] 



1 + . . . + k^s"^ + . . . ± 5^(1+^)/^ 1 + U^'^ + . . . ± 65^'^/^ ± • • 



±B'ts [I + 6;/^/^ + . . . ± y^s'''^ ± . . .] , (3.12) 

where, again, all the coefficients are noncritical and the same remarks as before apply to the 
signs of yl, X, and B. The correction factors for the B^ B and B' terms run through 

all powers of s with exponents of the form (no + 9[h])/ (3 (recalling the definitions ( |2.16|) , 



( p.24| ), etc.); terms with odd hi carry ± signs; when no = we have a/^ = a/^^ and b/^ = b/^. 
Expressions for ^4, A, etc. are given in (|4.28|) - (|4.31|) . 



Now note that the amplitude B^ carries a factor t which vanishes at Tg. Away from the 
endpoint this term induces a linear variation of rh with m which simply means that the 
tangents to the binodals at the A point (for T ^ T^) are no longer parallel to the tangent 
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at the endpoint. Such a variation is, of course, to be expected and does not represent any 
real change of shape as T deviates from Tg. To see this more exphcitly, note that we may 
redefine the coefficient Lp, which enters the definition (|3.3|) of m, as a noncritical function, 
Lp(t), chosen so that the tangent at the A point is always parallel to m = (i.e., to the m 
axis); then one has B' = Q while the other terms in (|3.12| ) do not change form. With this 
understanding for t 7^ we may conveniently define 

Am = m — mx{t) and Am = m — rhx{t)^ (3.13) 

which reduce to m and m, respectively, at the endpoint. 

The diameters of the nonsymmetric A binodals may now be found parametrically by 
multiplying out in (|3.11|) and ( |3.12|) and dropping all terms which carry ± signs. If the 



parameter s is eliminated in favor of x = Arhdiam/^i the diameters can be written 

Arridiam = Ax[l + /Cax^/^^"") + • • • + a.x'^/^'-''^ + • • •] 

+Wax(^+''s)/(^-") [1 + ■ ■ •] , (3.14) 

where we suppose a > while 

]Cx = {AK - Ak)/AA and Ux = Bh^. (3.15) 

We see that the slope {dm/ dm) of the diameter remains finite at the endpoint but, in 
general, the curvature diverges at the endpoint. 

The slopes = Am/ Am of the tielines follow similarly from the terms in ( p.llj ) and 
( P.12| ) carrying the ± signs. Using, again, x = Afhdiam/A as variable one finds, for a > 0, 

K 



_ ~a/(l-a) ^ _ ~^4/(l-a) ^ . . . 



{l-a)A 

^ (3.16) 



B 



As was anticipated, the tielines do not, in general, remain parallel to the A-point binodal 
tangent; however, the variation in slope is evidently slower than linear in Am. 
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Finally, one may eliminate s between and ( p.l2| ) directly and write the general, 

nonsymmetric A hinodals in terms of x = |Am/S| as 

+kx^l^ [1 + h.x'^l^ ± IaAx^''-^^'^ ± IkKx^^-^^'^ + ■•■], (3.17) 
where the ± signs refer to Am ^ (for B > Q) while 

aA = ax = -(1 - oi)/f3B, l/A, 

a4 = 54 — (1 — a)h4^/ (3^ 64 = —64//?, • • • . (3.18) 

We see that the leading behavior of the binodals, with exponent {1 — a)/ (3 (for a > 0), is 
the same as in the symmetric case (|3.1U| ). The surprising new feature, however, is the large 



number of numerically similar low-order correction terms. If we write the expansion for a 
general binodal in the form 

Am = ^^f|m|^% (3.19) 

i 

(with ± for m ^ 0) the nonsymmetric A-hne binodals generate the exponent sequence 

^1(5 = 1 - tt, 1, 1 + 1 - tt + ^4, 2 -2a- (3, 1 + ^4, 2 - « - 

2-/3, + •••• (3.20) 

For (i = 3 the Ising numerical values are 

i;ll3 ~ 0.895, 1, 1-328, 1.44, 1.46, 1.54, 1.57, 

1.67, 1.9, (3.21) 

where, here and below, we use the rough approximation 9^ 1.0; when ^ 4— one gets 1, 
1, |, 1, |, 1, |, |, |, Finally, note that the presence of the various ± signs in ( |3.17|) 
reflects the nontrivial behavior of the diameter and consequent lack of binodal symmetry 
outside the innermost asymptotic region. 
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C. Spectator Phase Boundary at the Endpoint 



The spectator phase, a, is bounded in the space of thermodynamic fields by the surface 
a (see Figs. |l]and H) which may be specified by the function ga{t^ h) which, as explained in I, 
is found by equating the a and /37 free energies, i.e., by solving G^{g(j^ t, h) = G^^{ga^ t, h). 
In leading order this was carried out in I but, for the present purposes it is useful to have 
the results correct to higher order. Here we present expressions for T = (or t = 0), i.e. 
on the endpoint isotherm. 

The detailed analysis is presented in Sec. |V|.C where one sees that it is advantageous to 
retain /i as a principal variable. The results for the symmetric case are the simplest in form: 
we find 

g^it = 0;h) = -J\h\^'+'^/'Zs{\h\) - J2h' - Mh\'^^"'^ + • • • , (3.22) 
where the singular correction factor is 

Zs{z) = 1 ± Ciz(^-")/^ + C2^2(^-«)/^ + C3^2-(^/^) 

+ c,z'*/^ ± c'^z^^--+'*^/^ + c,z'+^''/^^ + ■■■. (3.23) 

The upper (plus) signs in Zs correspond to case B or gi < 0; recall ( |1.4|) and Fig. ^; the 
lower (minus) signs describe case A when qi > 0: see ([T]^) and Fig. [T|. 
The leading amplitude in (|3.22| ) is given, using (|2.26| ), by 



J = QM(^+^y^W^/(D^ - roD,), (3.24) 



where Qe and U are defined via ( |2.13|) and (|2.10|) while, for the symmetric case, one has 



roi^3 = and J > 0. In addition we state 

ci = J/t/'/^, 02 = wlqlJ^/U^'^, (3.25) 

while the other coefficients are recorded in Sec. |V|.C. The result ( p.22| ) can be expressed in 
terms of h by using 
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h = h[l- n T rici + . . .] (3.26) 

which, however, is vahd only for t = and g = g^. We note that {5+l)/5 = {2 — a)/ A 1.21 
in agreement with I: see also Fig. | for a portrayal of ga{0^ h). We defer discussion of the 
correction exponents until the noncritical/spectator binodals are presented: see Eqs. 
and (CT ). 



In the nonsymmetric case the leading variation of gcr{t = 0) is, in general, linear in h 
(and h): see Fig. H. Specifically, subject to 

J, = Ds/{D^-roDs)^0, oc, (3.27) 

we find 

g^{t = 0]h) = -Jji - J\h\^^^^^'^ZN{\h\) - J2h^ 

-Jji\h\^'^^)l' - J^\h\^^'^^y' + . . . , (3.28) 

where the norjsymmetric singular factor has the expansion 

+ + )^2(l-a)/A ^ ^^4^3-(3/A) _ ^ ^^^.^^2-(l+a)/A ^ rf;|'^2-(l/A) 

- a.id'l + d\d,)z^'^+'-"^'^ (3.29) 

The two signum factors are given by 

at = sgn(t) = sgn(g/i), ah = sgn(/i) = sgn(ji/i), (3.30) 

in which we suppose the coefficients 

q = qo- qi{Ds/D,), Ji = 1 - r^{D^/Dr), (3.31) 

are nonvanishing; this will be true in the general nonsymmetric case. (We do not analyze 
the exceptions although no problems of principle arise.) 
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We see from ( |3.28|) - (|3.30|) that terms which change sign are not now determined simply 



by the slope of the A line (case A or case B), as in the symmetric situation, but rather 
by more complicated considerations. This arises simply because the manifold t = in the 
(g^t^h) space, see Fig. |1|, can cut the plane t = in various ways. For small asymmetry, ji 
remains positive giving ah = sgn(/i) but q may be of either sign. As expected from I, the 
leading singularity in is the same as in the symmetric situation; however, the corrections 
now contain further, new powers. 

The leading correction amplitudes in are 

(ii = d[ = - qoro)J/U'/^. (3.32) 



The remaining leading coefficients are listed in Sec. |I^.C. As before the result (|3.28|) can be 
expressed in terms of h by making the substitution 

h = nh - j\h\('^'^/' + fh\h\'/' - a,f\hf-'<^-^y^ - j,h' + ..., (3.33) 

where j = roJjiljil^^^^"^^^ while /, etc., are given below in (|4.48| ). 



D. Noncritical Endpoint Binodals 

We are now in a position to answer Widom's question regarding the shape of the non- 
critical or spectator-phase binodals, i^g^, at the endpoint. The essential point is that the 
densities pi and p2 and, hence, m and m, are noncritical functions of g^ t, and h in the 
spectator-phase a since G^{g^t^h) is noncritical. Consequently, on the endpoint isotherm, 
t = 0, the singular shape of the a binodals directly reflects the singular shape of the phase 
boundary 5'cr(0, h). 

To state the results for the symmetric case we introduce the endpoint susceptibilities 

Xe = -2G9 > and Xe = -2G^ > 0, (3.34) 
and the endpoint density 
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= -C" + < 0. (3.35) 
The noncritical binodal is then given by 

m = < - - C2XI - Csa;^/' + ■ • ■ , (3.36) 



where 



X, = m/x:, C = JTe. C2 = D,re/Dr, C,= [^-%\ J^Xe , (3.37) 



while Zs{z) is given in (|3.23|) . The leading exponent is {8 +l)/8 = (2 — a)/ A as stated in 



the Introduction. If we use the general binodal expansion ( |3.19|) the sequence of exponents 
arising now is 

?/;fA = 2-a, + 3 - 2a, 4 - 2a - 2/?, 3-2^ + ^4, 4 - 3a, 

5 -3a -2/?, 4 -2a-/? + 05, (S), (3.38) 

with values for Ising = 3: 

V^f ^ I.2O9, 1.55, 1.78, 2, 2.13, 2.35, 2.57, • • • , 2.9, • • • , (S) (3.39) 

(using, again, 6^ 1.0). 

In the limit d —> 4,— the sequence for ^/^f is |, |, 2, 2, 2, |, |, • • •, |, • • •. Note that 
the leading correction exponent found by Klingeri was |. His classical phenomenological 
treatment should correspond to ^ 4— but | does not appear here: the reason is that 
he did not (expressly) consider the symmetric situation. We also find the exponent | (and 
others) when symmetry is lacking. 

In the nonsymmetric case, the endpoint susceptibilities become more complicated: we 
find they are given by 

x: = -2{G^ - 2L,Gt + LlG2) > 0, (3.40) 
x: = -2{Gt - 2T,Gt + r^G^) > 0, (3.41) 
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where the significance of the axis slope, L^, was explained in Sec. |T|.A above. From 
and (1^) we obtain 



L. = (G^-G^)/(G?-G?), 



(3.42) 



while the endpoint density is 



= ro(G^ - G°) - + < 



(3.43) 



Using, again, = /Xe a variable, the noncritical endpoint binodal in the non- 
symmetric case is expressed by 



(3.44) 



where ± corresponds to /i ^ while the amplitudes m", ?7i2 , • ■ • are presented below. The 
linear variation of with shows that the tangent to the noncritical endpoint binodal at 
the endpoint E'^ is, in general, not parallel to the tangent at (the m-axis): see Fig. |^(b). 
The corresponding amplitude, m", is 



< = 2(-G^ + roG^) + 2L,{Gl - roG^). 



(3.45) 



The leading singular exponent, namely 1 + (1/5), is evidently the same as the symmetric 
case, while the amplitude, m2 , is 



(3.46) 



Recall that the coefficients j, Ji, and J are defined above in Sec. |T|.C. 

The leading correction exponent is now just that found by Klingeri in his classical treat- 
ment; it does not appear in the symmetric case. The expression for its amplitude, , is 
complicated but, for the record, we quote the result, namely, 

+ /x? 



mg = 



+2sgn(ji)5f2 



<^-2(iJi- J)(-G2 + roG?) 



Tfl 



(3.47) 
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where the new coefficients are 



X? = 2{-Gt + L,Gt){jJi - J), (3.48) 
92 = rl\n\'^('/'^J\h - ^-^U^'jJ. (3.49) 



E. Critical Endpoint Binodals 

Now we conclude our discussion of the endpoint itself by presenting, finally, the shape of 
the critical phase binodals, i3f and B]. These can be obtained by using the thermodynamic 
potential for the critical phase, G^^{g, t^h)^ and the endpoint phase boundary, Qaih). Details 



are given in Sec. [I^.D. As discussed before, it is convenient to describe the binodals with 
the aid of a parameter s (in this case, related to 1^1^/"^) which vanishes at the endpoint and 
increases into the /? and 7 phases. 

In the symmetric case, S, the critical endpoint binodals may then be specified by 

m = ±Es[l + u^s^'/f^ + + • • •] ± Vis^/^[1 + vis^/f^ + • • •] 

^y^^il-a+AW [1 + . . .] ± y3^(2-«+A)//3 [l + . . .] ^ (3.50) 

m = ^^(^-"^/^[l + u^s^'/^ + uis^^-""^/^ + •••] + + ...], (3.51) 



where ± corresponds to /i ^ 0, and the coefficients are given in ( [4.58|) and (|4.59|) . 



The symmetric critical endpoint binodals may finally be expressed in terms of m as a 
variable by solving (|3.50|) for s and substituting in ( |3.51|) . With x = \m/E\ this yields 



m 



^^(1-^)/^ [1 + u^x^^/^ + uix^^-""^/^ + •••], (3.52) 



where 



Ui = Ui — {1 — a)ui/ jS^ = 'U4 — (1 — a)u4/p. (3.53) 

The term in E provides the dominant behavior with the same exponent as the A-line binodals 
given in Sec. |TTI|.B. One should notice that the amplitude E is negative in case A while it is 
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positive in case B owing to the negative sign of discussed following ( |2.23|) .ll3 Hence it has 
the same sign as the amplitude A of the lambda line binodals: see ( p.lO| ). This also holds 
in the nonsymmetric case. 

Indeed, the nonsymmetric^ N, critical endpoint binodals may be described similarly. In 
terms of the parameter s we find 

m = ±Es[l ± u,s^^-'^/^ + U2S^'-^^/^ + . . . + u,s'^/^ + • • • ± u,s'^/^ + • • •] 

± ^^,(A-i)//3 + ^2,(1-)/^ + . . . + ^^/4/^ + . . . ± ^^/5//3 + . . .] 

±V2S^/f'[l ±voS'"]+ Vss^'-^^/f'il + ...], (3.54) 
rh = Es^'--^/^[l ± u,s^^-'y^ + £^2^(1--)/^ + . . . + u,s'^/^ • • • ± u,s'^/^ + • • •] 

±Vis'^/^[l ±vos + ---] + y25^'""^/^[l + •••], (3.55) 



where the leading coefficients are presented in (|4.60|) - (|4.62|) . Solving for s in ( |3.54|) and 
substituting into (|3.55| ), one finally obtains 



rh = + u^x^'/^ ± ± 222x(^-")/^ + ...], (3.56) 
where the leading coefficients are 

ui = ui-{l-a){E/Vi + m)/f3, 

U2 = Vi/E, U4 = U4 — {1 — a)u4/P, (3.57) 

while the correction factor exponents have d = 3 Ising values 64/ f3 1.64, {A — l)/f3 
1.73 and (1 - a)//? ^ 2.73. 

F. Binodals above the Endpoint Temperature 

Let us consider, first, the spectator-phase binodal above Tg — see Figs. |^(c) and ^(c) 
— which is the simplest to analyze. Since G^{g^t^ h) is noncritical, the densities m and m 
are noncritical functions of t, and h in the spectator-phase a. At fixed t > 0, the phase 
boundary g^{t] h) is also a nonsingular function of h with t-dependent expansion coefficients, 
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which are discussed exphcitly below in Sec. |^.E. Consequently, on the isotherms above Tg, 
the a binodal becomes noncritical. However, singularities of the binodal are to be expected 
as T ^ Tq+. 

In the symmetric case^ by using the previous definitions (|3.34|) and (|3.35|) and the phase 
boundary ga{t] h) given below in ( |4.67| ), we obtain 



rh = fat- TeigU + - TeQU-^xl + • • • , (3.58) 

where = Tn/Xe^ as for the symmetric noncritical endpoint binodals in ( |3.36|) , while the 



coefficients, g^^^ etc., are given below in ( |4.68| ). Note that the curvature of the binodal 
diverges like when T ^ Te+. 

In the nonsymmetric case^ using ( [4.70|) , we obtain 

m = + m^x« + rh^{t)xl + • • • , (3.59) 



where rh^ and are given above in (|3.43| ) and ( |3.45|) , respectively, while the coefficient of 
second order in (= m/Xg) is 

= 2iG2 - L,Gt)9X^jlt-'' + ■ • ■ , (3.60) 

where g^^ is given below in ( |4.71| ). Here we have neglected higher order corrections in t. 
Just as in the symmetric case, the curvature of the binodal diverges when T —> Te+. 

Consider next the critical phase binodal B^^ above Tg: see Figs, ^(c) and ^(c). This may 
be determined using (|4.15|) below and its twin for m with the aid of the spectator-phase 



boundary, ga{t] h)^ which is derived in Sec. ^.E. For fixed t > 0, the small y expansion for 
the scaling function W^{y^y4^ • • •) yields only integer powers of h in (|4.15| ) and its twin so 



that the densities m and m are noncritical functions of h. Consequently, the critical phase 
binodal is again noncritical above Tg. 

In the symmetric case^ the densities m and m can be written in terms of h by using 
( |4.15| ) and its twin as 

m = ht-^h + "', (3.61) 
rh = mo{t) + kt-^-^h^ + • • • , (3.62) 
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where fh{){t) is a function of t only while the coefficients are 

h = 2Q,U^Wl^\l-q^{D2/D,)\-\ 

h = -jqiQeU'Wl,\l - qi{D,/D2)r-\ (3.63) 

Notice that I2 is negative in case A while it is positive in case B, as for ^4, the leading 
amplitude of the lambda line binodals: see the paragraph below ( ^^ . We may eliminate h 
between ( p.61|) and (|3.62| ) and write the binodal in terms of x = m//i, noticing li > 0, as 

m = mo{t) + ht^-^x^ + "'. (3.64) 

Since 7 > 1 in the d < 4 Ising universality classes, the coefficient of the quadratic term in x 
vanishes as T ^ Tq+. This result could be anticipated, since the critical endpoint binodals 
have the leading exponent (1 — a)/ jS [c^ 2.73) in the symmetric case. Thus the curvature of 
is singular but nondivergent when T ^ T^. 

In the general nonsymmetric case^ the situation is more complicated. The densities can 
now be expressed as 

m = mo{t) + hr^h + hr^-^h^ + • • • , (3.65) 
fh = mo{t) + ht^-'^h + y-^-^h^ + • • • , (3.66) 



where the constant coefficients are presented below in (|4.74|) . Note that the term linear in 
h for m has a leading t-dependent coefficient that vanishes when T —> T(.-\-. As before, the 
critical phase binodal can be written in terms of x = Am/li with Am = m — mo{t) as 

m = rho{t) + + l2t^-^x^ + • • • . (3.67) 

Evidently, both the coefficients of x and x^ are singular but vanish when T ^ Te+ and 
7 > 1. 



G. Binodals below the Endpoint Temperature 



Below the endpoint temperature three phases, a, /?, and 7 may coexist on the triple line 
r. The binodals near a triple point then spring from the corners of a three-phase triangle. 
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The corresponding phase diagrams in the density plane are shown in Figs, ^(a) and ^(a) for 
the two cases NA and SB, respectively. Thermodynamic stability then requires that these 
diagrams must satisfy Schreinemakers' rulesiiE^ details are given in Ref. 17. 

The explicit forms of the spectator-phase binodals^ ^<^, can be obtained without difficulty 
by using the phase boundary gcr{t^h) below Tg as presented in ( [4.76|) and (|4.79| ) for the 



symmetric and nonsymmetric cases, respectively. In the symmetric case^ the binodal may 
be expressed as 

where = rn/Xe the upper (lower) sign corresponds to m > (< 0), while the 
coefficients, g~Q^ etc., are given below in ( |4.771 ). Note that the slope vanishes as T ^ Tg— 
while the curvature diverges as \t\~^. In the nonsymmetric case^ the binodal is given by 

m = + m^{t)xa + rh^{t)xl + • • • , (3.69) 

where the coefficients are 

= < T 2^-2Ji(-G'^ + L,Gt)\tf + • • ■ , 

m^it) = -2g-^jl{-G1 + L,G'«)|t|-'^ + • ■ • , (3.70) 

while the =F signs again correspond to m ^ 0. The coefficients, g~^Q^ etc., are given below 
in (|4.80|) . Notice that the linear terms do not vanish, but approach the same value when 



The critical phase binodals^ and can be obtained, in principle, by using ( [4.15|) 



and its twin with the aid of the phase boundary g^{t^ h) given below in ([4.76|) and ( ^.79|) . 



However, the analysis becomes more complicated, since these binodals are associated with 
the lambda line binodals near the vertices of the three-phase triangle: see Figs, ^(a) and|^(a). 
Hence, we do not present their explicit forms here. One can anticipate, however, that the 
binodals have linear slopes and quadratic terms which both vanish when T ^ Tg— in the 
< 4 Ising universality classes. 
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IV. DERIVATION OF THE BINODAL EXPRESSIONS 



In this section we sketch, for completeness, some of the details entering the derivation of 
the results for the binodals presented in Sec. |TTl| from the postulates of Sec. In addition, 
we give explicit expressions for the leading amplitudes entering the formulae of Sec. |T| in 
terms of the original parameters of the postulated free energies of Sec. |I|. 

A. Principles for Obtaining Isothermal Sections 

Our aim is to describe isothermal sections of the full {g^ t, h) phase space in terms of the 
density variables 

pi = -dhG, p2 -dgG with dh = d/dh, dg = d/dg. (4.1) 

Accordingly, we treat t as a fixed parameter and regard only g and h as varying. The basic 
nonlinear scaling fields i and h are then to be viewed as functions only of g and h. Once the 
appropriate derivatives with respect to g and h have been performed, however, it is more 
convenient, in light of the scaling postulate ( |2.13|) , to employ the nonlinear scaling fields 



i and h as the primary field variables. Note, in particular, that both the A line and the 
triple line, r, lie in the plane h = 0. Beyond that, the A-line or p-surface binodals also 
correspond to h = while the spectator-phase and a binodals are of interest only for small 
h. Consequently we express g and h in terms of i and h via the noncritical expansions 

9 = 9x{t) + eii + e2h + e^P + ejh + efh^ + • • • , (4.2) 
h = hx{t) + ht + hh + /3P + hih + hh^ + ---, (4.3) 

where the A-line values, gx and /ia, were introduced in (|2.8| ) and seen to be noncritical 
functions. Likewise, all the coefficients ej{t) and fj{t) are noncritical with, in the symmetric 
case, 

S : 62 = 64 = A = /a = /s = 0, 

ei = + 0(t), h = 1, 63 = 65 = h = (4.4) 

qf qi qi 

29 



More generally, with Aq = qi — ro^o (7^ 0), we have 

N : ei, 62, A, /2 = (1, -go, -ro, gi)/Ao + (4.5) 

while 63, • • • , /s are also readily found in terms of the qj and rj. 

Any noncritical property P{g^t^ h) with expansion can then be rewritten as 



P{g, t, h) = Px{t) +Pi{t)i+P2{t)h +Ps{t)P + . . . , (4.6) 
where the value on the A line is given by 

P^{t) = P, + Pxit + Px2t' + ---^ (4.7) 

PXI = PlAg + P2 + Ps^h. PX2 = PlAg2 + ' ' ' + P^ ' ' ' , (4.8) 

where A^, A/^, A^2, etc. are defined in ( |2.8|) and (|Z9|) , while the remaining noncritical coeffi- 
cients take the form 

Pj{t) =Pj, +Pjit +Pj2t^ + . . . , = Piej + Psfj. (4.9) 

Pji = 2(P4A,e, + P.Agfj + P.AhCj + Pee, + + i^gA/,/,), (4.10) 

for j = 1 or 2, and 

Pse = Pies + P3/3 + P^ej + 2P5ei/i + P9/1 , (4.11) 

and so on. 

Of course, we eventually wish to eliminate i and h in favor of pi and p2 or, in view of 
the discussion of Sec. |ITI|.A in terms of 

m = -dG + {dG)e and m = -dG + {dG)e, (4.12) 

where the compound differential operators are 

d = dh-L^dg and d = dg - L^dh. (4.13) 

However, once we have expressions for m and m in terms of i and /i we can regard these fields 
merely as auxiliary parameters relating m and m. Note in particular, that coexisting phases 

30 



must have the same values of i and h. Thus for the p binodals we can put s = (— t)^, for 
t < 0, and set h = 0. This indicates the origin of the parametric descriptions of the binodals 
presented in Sec. |Tl|.B. Similarly, for the binodals associated with the a phase boundary, 
equating the free energies G^^ and gives a relation for i in terms of h (and t); then h is 
an appropriate parameter. 

The axis slopes L^r and Lp in ( |4.13|) were explained in Sec. |ITI|.A and the slope L^r was 
given in ( |3.42|) . Below we will establish the t-dependent result 



Lp{t) = ro + [2nAg + rs + viAh - r^inkg + T2 + 2r^Kh)]t + 0{t^), 



(4.14) 



where Lp{t) was introduced just before ( p.l3|) and Lp = Lp{0) = tq. 

Now using (|4.12| ) and (|2.13|) we obtain the primary density in the form 



m = {dG^)e - dC^ + \i 



2-a 



{dQ)W± + J2{dUk)QWi'^\if^ 



k>4 



±{di)\i\^-''QW± T {di)h\if-^AQUW^ + {dh)\ifQUW^ 



±5 



(4.15) 



for /i ^ 0, with a precisely similar expression for m with d replacing 5, while 
Wi{y, 1/4, ■■■) = i2-a)W± + J2 ^kUk 



k>4 



{dW±/dy), Wi''\y, ■ • ■) = {dW±/dyk). 



(4.16) 
(4.17) 



Note that dG^ and the coefficients dQ^ dUk-> di^ and dh are all noncritical and so can be 
written as in ( [4.6|) . This form thus enables one to identify all the singular terms appearing 
in m and rh. 



Now on the A line we have t = h = 0. Thus (|4.15|) and its twin for rh yield the expansions 
and (|3^ ) for rux and rhx with 



Ml = 2 [L,(A,G° + AhGl + G°) - AgG^ - - A^G^] , 
Ml = 2 [L.iAgGl + G? + A;,G°) - A,G° - A^G^ - G°] , 



(4.18) 
(4.19) 



SO that Ml = and Mi = 2(^4/^1 — Gq) in the symmetric case. Defining R{g^t^h) 
{dG^)e — dG^ and R likewise, and expanding as in (|46| ) yields, for j = 1, 2, 
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= 2 [LAGh + Glf,) - G°e, - G^fj] + 0{t), (4.20) 
4 = 2 [ro(G°e, + - G°e, - G°/,] + 0(t), (4.21) 



where (|4.14| ) was used for Lp. For reference below we also record 



{di)x = qo- L^qi + feA^ + ^qe^h + qi - L^{2q2Kg + q^ + q^Kh)] t + • • • , (4.22) 

{dh)x = 1 - L^ro + [riAg + rs + 2r3A/, - L,{riAh + 2r4A^ + r^)] t + • • • , (4.23) 

(at)A = gi - mo + [2g2A^ + 93 + gsA/, - Lp{q^Ag + 2g6A7 + qr)] t + • • • , (4.24) 

(a^)A = ro - + [nAh + 2r4A^ + rs - ^^(riA^ + + 2r3A/,)] t + • • • . (4.25) 

Clearly, any desired higher order terms in the t, h expansions can be obtained straightfor- 
wardly. Finally, we remark that we will shortly see that the condition determining Lp{t) is 
that {dh)x vanishes term by term; substitution of ( |4.14|) in (|4.25| ) checks this. 



B. Derivation of the A-Line Binodals 



The binodals associated with the A line may, essentially, be obtained directly from ( [4.15|) 



and its twin, by letting h 0± with i < 0. In doing this the small y expansions (|2.21| ) must 
be used with attention to the (Jhi{y) factors defined in (|2.22| ). When this is done the 
terms in (|4.15| ) generate only integral powers of the terms in act merely to modify 



the correction factor of the term. Note that the i and h expansions of Q and of the Uk 

yield correction terms varying as Itl^+^I'^l for all integers n > and all > 0. The term in 
h\i\^~^^ which diverges as |t| ^ 0, vanishes identically. Lastly, the term in contributes 
both to m and rh. 

Introducing the parameter s = then yields the previously quoted expansions (|3^ ) and 
( p.7|) for m and m in the symmetric case. The linear term in s is absent in this m expansion 
because the coefficient (dgh) vanishes identically by symmetry when h = h ^ and Lp = 
is dictated. Similarly, terms varying as s^^^ and s^^~^'^/^ are absent in the expression for m 
since L^- = and thence dG^^ dQ and di all vanish. For even k the derivatives dhUk = dUk 
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(for L(j = h = 0) also vanish by symmetry. However, in the fuUy symmetric situation each 
odd scahng field, U2j+i{g^t^h)^ must itself be odd in h: see (|2.12| ). Hence after operating 



with dh^ contributions with odd k in the terms proportional to |tp ^^^^{dUk) in ( [4.15|) 



appear in the expansion for m in the symmetric case. Since 2 — a = /? + A these terms are 
responsible for the appearance of the correction factors = m (|3.6|); see also (|3.9|). For 
completeness we record the leading amplitude values 

A, = -{2-a)qrQM\. B, = UQM-i. Ke^^Gl/qr, (4.26) 

Clearly all other amplitudes are readily generated although their complexity increases rapidly 
with order. 

In the general nonsymmetric case the Uk for odd k need not vanish when h ^ but 
the scaling function, W-{y^ ?/4, 7/5, • • •), still has special behavior for small yk when k is odd: 
see (|2.18| ). This is the reason why the ± signs (corresponding to h ^ 0±) appear in the 



expansion ( |3.11|) for m. The expansion for the secondary density m, when initially generated, 
has a similar structure. In particular, the leading term is proportional to = s. However, 
at this point we should, as explained in Sec. |Tn|-A, complete the specification of the density 
m by appropriately choosing Lp{t). This should be done by examining the common tangent 
to the critical binodals, namely i3f and B^^ at the endpoint: see Figs. H(b) and ||(b). But 
these binodals involve the a phase boundary which we have not yet studied. Instead, we 
will select Lp so that the common tangent of the A-line binodals B^^ and B^~ or B^ and 
B^ coincides with the m = axis when extrapolated to the endpoint. It will be confirmed 
below that this criterion gives the same value for Lp. The coefficient of the offending 
term is {dh)^'. see ( [4.25|) . This vanishes when Lp = tq so confirming (|4.14|) for t = 0. 



The residual t and \i\ dependence of (dh) then yield the B^ts and Bs^'^^^'^^^ terms in the 
expansion ( |3.12|) for rh. The latter term is unavoidable in general and further complicates the 



singular corrections to the p binodals in the nonsymmetric case. Nevertheless, as explained 
in Sec. |lTl|.B, the former term, linear in 5, can be eliminated by adopting a temperature- 
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dependent definition for rh by allowing Lp to vary noncritically with T. The criterion now 
is to make {dh)x{t) vanish. Reference to (|4.25|) then confirms the leading term in Lp{t) 



presented in ( ^.14|) . 

The leading amplitudes in (|3.11| ) and ( p.l2| ) for the nonsymmetric case are, recalling 
and ( ra , 

A{t), A{t) = -(2 - a)QxW%[{di)x^ {di)xi (4.28) 

B{t) = QxUW'_,{dh)x^ (4.29) 

= -Q^UW\[{2n - rori)ei + (n - 2rors)fii (4.30) 

Ke = -i?le, Ke = -Rle, ^4 = 64 64. (4.31) 

One further has = a^^ = ^5, etc., although correction terms carrying 'noncritical factors' 
s^/^ = \i\ do not in general satisfy corresponding equalities. 

C. Spectator Phase Boundary: Endpoint Isotherm 

As indicated in Sec. fTJ-C, the first step in studying the binodals not associated with the 
A line is to obtain the phase boundary a as specified by gcr{t^ h). On recalling (|2.26| ) and 
( p.l3| ), one sees this is to be found by solving 

D{g,t, h) = -Q\i\^-''W±{y,y,, ■ ■ •), (4.32) 

where D{g^ t, h) is noncritical with 0^ = and Di > 0. Here we focus only on the endpoint 
isotherm, T = Tg or t = 0. Now consider the argument y in leading order, using ( |2.5|) and 

y = Uh/\i\^ ^ U{h + rog)/\qoh + qig\^. (4.33) 

If To, Qo and qi do not vanish (as in the generic norjsymmetric case) it is evident that when 
g, h ^ on a one in general has y ~ [max(|5f|, l/i])]^""^ which diverges to oo since A > 1. 
Thus to study ( |4.32| ) on the endpoint isotherm we must utihze the large y expansions ( ^.231 ) 
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for the scaling functions entering (|2.15| ). In the symmetric case one actuaUy has Vq = = 0; 



but it then transpires, as shown below, that ^ \h\^'^ ^^^^ so that y ^ \h\^ ^ . Since a <1 
we see that y again diverges. Thus in (|4.32|) we must always use the expansion 

+W^^y,\y\^'-'^^''^/^{i ± wf\y\-'/'' + ■ • ■) 

+W^k.sgn{y)\y\^'--+'^^y^{l ± wf^\y\-'/^ + • • •) 

+ •••, (4.34) 

where the ± signs correspond to t ^ 0. 

The analysis is considerably simpler if one uses /i as a variable in place of h. To this end 
we rearrange and with t = to obtain 



h = h - r^g - {ri - 2rQr^)gh - rji^ - r^g'^ H , (4.35) 

where = — TqTi + r^rl^ and 

i= qoh + Pig + p2h^ + Psgh + p^g'^ H , (4.36) 

where the leading coefficients are 

Pi = qi- go^o, P2 = q6- qors, 

P3 = q5- qon + 2qorors - 2qQro, Pa = q2 - q^ro + ge^o " 90^4- (4.37) 

Note that in the symmetric case one has qo = q^ = qr = 0^ Tq = = = and so ps = 0] 
we may suppose 7^ 0. 

Now, combining these results for the symmetric case yields the asymptotic equation 

D,g = -D,g' - D^h^ - {Q, + Qig + • • W^'^'^^'Z - • • • (4.38) 

with scaling factor, from ([4.34|) , 



+W^^U,ig, 0, h)\Uhf^/^[l ± w^^ \y\-'/^ + • ■ •] 

+sgn{y)W^')U,{g, 0, h)\Uh\''/^[l ± ■•■] + ■•■ . (4.39) 

35 



These equations are to be solved together with 

1^1 \(li9\ 



\y\ 



-l/A _ 



l + —g + 

Qi Qi g 



(4.40) 



\uh\y^ \uh\y^ 

to yield g = ga{h). This can be accomplished iteratively by noting that in leading order 
g^ fa —J\h\^^^^'>/^ ^ where J was defined in (|3.24| ); however, care is called for! 

One obtains the result quoted in (p.22| )- (|3.25| ) which may be supplemented by 



J2 = D^/D^, J, = [{D,/Di) - {Qi/Qe)]J\ 



(4.41) 
(4.42) 
(4.43) 



where U^^s is the first nonzero expansion coefficient of U^^g, 0, h) ~ U^^^h in the symmetric 
case. The expression (|3.26| ) for h{h) on a follows from (p.6| ) and ( p.22| ) by reversion. 

The phase boundary in the nonsymmetric case follows in the analogous way but greater 
care is needed because of the increased number of nonvanishing and competing terms. Thus 
on the right side of ( ^4.38|) the new terms —D^h and —2Dr,gh appear, where Dr, = Dr, — 
\Dz{ri — 2rQrz) — DqTq. The former term dominates and so in leading order one now finds 



ga ~ -Jih - J\h\ 



{5+l)/S 



(4.44) 



where Ji was defined in (|3.27|) . This in turn yields the new behavior 



\y\ 



-l/A 



k1 



1^1 



1-(1/A) 



+ 



(4.45) 



q q 

where q = qo — piJi = q/ji while q, ji, and dh were defined in ( p.30| ) and (|3.31|) . 

In this way one obtains the result (|3.28|) - (|3.32|) which must be supplemented by new 
expressions for J2 and J3 while 



(4.46) 



The expressions for dg, dg, are long and uninformative but we quote 
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4 = wf^d5\q\/U^/^, d'l = wf^d^p^J/U^'^. (4.47) 
Finally the remaining coefficients in (|3.33|) are 



32 = ji[roJ2 + (ri - 2rors)Ji - - f^Jl]. (4.48) 



D. Derivation of the Critical Endpoint Binodals 

The critical phase binodals at the endpoint may be obtained from (|4.15|) and its twin 
using the endpoint isotherm, ga{h)^ obtained in the previous section. In order to do so, it is 
more convenient to rewrite ( |4.15D as 



\2-a 



{dQ)Wi + J2idUk)QW'f^\if- 

fc>4 



±{di)Q\i\^-''W± + {dh)UQ\ifW^^, (4.49) 

and similarly for m with d replacing 5, while 

W±=W±- AyW^, (4.50) 

where A = 2 — a — P has been used. At the critical endpoint, t = 0, we use h as an 
auxiliary parameter relating m and rh. Using (|4.35|) and (|4.36|) , the noncritical functions, 
(9G^), (5(5), etc. can be expressed in terms of h. Recalling the general expansion for 
a noncritical function P{g^ t, /i), we find, for t = 0, 

P{g, t = 0, /i) = Pe + Psh + (Pi - roPs)ga + • • • , (4.51) 
and similarly for the derivatives 

dP = Ps- L^Pi + 2{Pg - L„P,)h + 2(P5 - L^Pi - roPg + roL,Ps)ga + ■■■, (4.52) 
dP = Pi- LpP3 + 2(P5 - LpPg)h + 2(P4 - LpPr, - roPs + roLpPg)g^ + ■■■. (4.53) 
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Likewise, in terms of gcr{t = 0, h) we obtain 

di= qo- L^qi + {2qQ - L^q^)h + {q^ - 2L^q2 - 2roqQ + L^rQq^)g^ H , (4.54) 

dh = l- L^ro + (2r3 - L^ri)h + (ri - 2L^r4 - 2r^r^ + L^r^ri)g^ H , (4.55) 

di= qi- LpqQ + (gs - 2Lpqe)h + (2(?2 - L^q^ - r^q^ + 2rQLpqQ)g^ H , (4.56) 

dh = ro - Lp + (ri - 2Lpr^)h + (2r4 - L^ri - rori + 2roLprs)g^ H . (4.57) 

As discussed before, the argument y of the scahng functions W± diverges to oc when the 
endpoint is approached on the a surface. Thus in ( [4.49| ) and its twin the large y expansions 
( p. 23 ) for the scahng functions must be used with attention to the (7/^{y) factors defined 
in (p.22|) and the multiexponents 9[k] in ( p.24|) . When this is done, we finaUy obtain the 
critical endpoint binodals from ( |4.49| ). 



Introducing the parameter s = {hl^^^^ then yields the previously quoted expansions ( |3.50|) 



and (|3.51|) for m and rh in the symmetric case. The linear term in s is absent in the expression 



for m when we choose l/p(0) = Tq which reinforces previous results. In the expression for 
m the term in (|4.49| ) provides a linear term in h that yields the s^^^ term in ( |3.50|) ; 
the terms in provide the 5(2-«+a)//3 ^qyir and higher order corrections, since (dQ) 

generates h in leading order; the term in provides the 5(1-^+^)/^ term for (di) for the 

same reason; then, finally, the term in provides the leading s behavior. In the expression 
for m all the terms, except for one in 1^1"^"^, provide correction terms, 5^^"^^/^, in ( |3.51|) ; the 
leading behavior, 5(1-^)/^, is generated by the term in 
The leading amplitudes are 

E=[{2- a)/ A] E = qiQeW^w^.U^'-^^/^. (4.58) 

For the record, we also quote 

(2 - a) 

u, = -^j-^y,q.JlU^'\ = ^^^U,M'^/\ (4.59) 

yZ CM.) VV^W-^ 
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^1 = -2^^^' ^ = 2G^^ + [Qi + (2 - a)rige/A]H/^f/(^-'^)/^. 

In the general nonsymmetric case, the hnear term in s is stiU absent in the expression for 
m: see (|3.55| ). The expression for m in terms of s is given in (|3.54j) ; the terms in ( |4.49|) 



yield the s^^^ term, as in the symmetric case, while the terms in |tp ^ provide the s^^ ^^/^ 
term and that in \t\^~^ gives ^(i-^)/^; the leading term, 5, is still provided by the term in 
\t\^. In the expression for fa the leading behavior is s^^ ^ as in the symmetric case, again 
provided by the term; the term yields corrections of leading order 5^/^, while the 

terms in and give the s^^~^^l^ term in ( p.54| ). The required amplitudes are 

E = [(2 - «)/A] (1 - L.ro)geH^^f/('-'^^/^, 

^ = - ro(Zo)QeW^^ti^?f/^'-"^/^. (4.60) 
For the record, we also quote the correction amplitudes 

= (go - L,gi)geVF^«^?f/('-"^/^, 
^2 = -2(G° - + 2Ji(G° - L,G\ - r^G% + roL,G°), 

1^3 = (Qs - L.gi)W^^/7(2-")/^, (4.61) 

= -2{G% - r^G%) + 2Ji(G'° - 2roG'° + r^G^), 
V^2 = [Qi + (2 - a)(ri - 2ror3)ge/A]l^^f/(2-")/^, 

and the leading further coefficients 



E. Spectator Phase Boundary: Isotherms above Tg 

In Sec. |V|.C, we studied the endpoint isothermal phase boundary, ga{h)^ in order to 
discuss the endpoint binodals. By the same token we study the phase boundary g^{t^h) 
above Tg as the first step in determining the supercritical binodals. This boundary is found 
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by equating the free energies, G^{g^ t, h) and G^^{g^ t, /i), of the spectator and critical phases, 
respectively, which yields ( |4.32|) with t > 0. The extended triple line f — see Figs. |1] and || 



— is defined by /i = for t > 0, implying ^ = 0. Since we consider only the vicinity 
of the extended triple line f, we must utilize the small y expansion (|2.19|) for the scaling 
function VI/+(?/, ?/4, ys, • • •) in (|4.32|). Using h as the principle variable, which is advantageous 



in discussing the critical phase binodal, the scaling function VI/+(^, ^4, 7/5, • • •) can be 
expanded in integral powers of h with t-dependent coefficients. The noncritical function 
D[g^t^h) can be expanded similarly. Then, solving ( |4.32| ) for g^{t]h) yields the desired 
nonsingular expansion. In this section, we consider only the leading t-dependent behavior 
of the resulting coefficients. 

Accordingly, we rearrange ( p.5|) and (|2.6| ) for t > using just the linear terms to obtain 



h = h- r_it - ro^ + • • • , (4.63) 
t = (1 - qor_i)t + q^h + {qi - qorQ)g H . (4.64) 



The higher order terms in (|2]5|) and ( ^^61) enter only as correction terms in the t-dependent 
coefficients. The noncritical function D{g^t^h) is then expanded, by recalling ( ^^21) and 
0, as 

D{g, t, h) = {Dr - r^D^)g + {D^ - r_rD^)t + dJi + ---. (4.65) 

Now we are in a position to find the isothermal boundary g^it] h) above Tg. 
In the symmetric case^ we obtain 

Drg + D2t + '" = -QW%\i\^-'' - QW%U^\i\-^h^ + 0(^^). (4.66) 

By symmetry only even powers of h appear. Solving for g with the aid of ( ^i.64|) then yields 



g,{t- h) = -gy - gy'-" - gyh' + • • • , (4.67) 

where the coefficients are 



7+3 = QDr'Wl,U'\D, - q^D^r. (4.68) 
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Notice that the coefficient of the quadratic term in h diverges as T ^ Te+. In terms of 
/i, which is advantageous in deriving the spectator-phase binodal, we obtain the same 
leading t-dependent coefficients for ga{t] h). 

In the nonsymmetric case^ terms linear in h appear in the expansion of the scaling 
function W^{y^ 2/4, 2/5, • • •) owing to the odd k exponents in ( |2.19|) . However, these terms only 
provide correction terms to the leading t-dependent behavior. Combining all the previous 
results yields the equation 

= -gi^+ol^T"" - QWl^U^\i\-^h^ + . . . . (4.69) 

Solving for g yields 

g,{t- h) = -g+,t - g+,e-- - jji - g^^.t^'h' + ■■■, (4.70) 
where Ji is given above in ( p.27|) while the other coefficients are 



g+o = {D2-r.iDs)/{D,-roD,), 

in which the numerical factor is 

= (1 - gor-i) - (gi - qoro)[{D2 - r_,Ds)/{D, - roDs)]. (4.72) 

Notice, again, that the coefficient of the quadratic term in h diverges when T Te+. The 
result (|4.70|) can be expanded in terms of h by making the substitution 



h = jih - r_it + • • • , (4.73) 



where ji is given in (|3.31| ). By utilizing (|4.70|) , the coefficients /i, • • •, /2 in ( |3.65|) and ( |3.66|) 
are found to be 

h = 2{l-L,ro)QeU^Wl^m-\ 
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h = -7('?o - L^qi)QeU'W%\lr-\ 

h = (2 - a)W%m^-''[{{q5 - 2rog6) - 2Ji(g2 - + r^qeWe 

+ {Qi - roqo){Q3 - MQi - r^Qs)}], (4.74) 
k = -liqi - roqo)QeU'W%\l\-'^-\ 

where Ji and to- are defined in (|3.27| ) and (|4.72|) , respectively. 



F. Spectator Phase Boundary: Isotherms below Tg 

The spectator-phase boundary, ga{t^ h)^ below the endpoint temperature can be obtained 
as in the previous section by using the expansion (|2.21|) for the scaling function W-{y^ 1/4^ • • •) 
in ( |4.32|) . The \y\ factors in ( |2.21|) yield the two branches of the phase boundary, gcr{t^h)\ 



see Figs. |^(a) and §(a). 

In the symmetric case^ combining the results in the previous section with the expansion 
yields 

Dig + D2t + -'- = -gi^^ol^T"" - QW^^U\i\^\h\ - QW%U^\i\-^h^ + • • • . (4.75) 



Solving this for g with the aid of (|4.64|) provides the result 

g^it; h) = -g-, - ^"iltp"" T g-2\tfh - g-,\trh' + ■■■, (4.76) 
where the upper (lower) sign corresponds to /i > (< 0) while the coefficients are 



g-, = D2/D,, g-, = QD^'+^W%\D, - qiD^]'-'', 

g-, = QD-^-^W\U\D, - q,D2f, g-, = QDl-^W%U^\D, - q.D^]-^ (4.77) 

Notice that the hnear term in h vanishes as T ^ Tg — , while the coefficient of the term 
diverges. In terms of h we obtain the same leading t-dependent coefficients for ga{t] h). 
Finally, in the nonsymmetric case we obtain the equation 

(Di - roD:,)g + {D^ - r_^D:,)t + D^h + ■ ■ ■ 

= -QW%\i\^-'' - QW^^U\if\h\ - QW%U^\i\-^h^ + • • • . (4.78) 
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By using ( |4.64| ), we can solve this for g to obtain 

g^{t- h) = -g-,t - g-,\t\'-" - (Ji ± 9^^)'^ - g-M'^h' + ■■■, (4-79) 
where, again, the upper (lower) sign corresponds to /i > (< 0), while the coefficients are 



(Di-roDs)' ' ^^'^ (Di-roDs)' 
Notice that the linear term in h does not vanish in the nonsymmetric case^ but the slopes of 
two branches approach the same value as T ^ Tg— . The coefficient of the quadratic term 
in h diverges as the endpoint temperature is approached from below. As before the result 
( |4.79| ) can be expressed in terms of h by using (|4.73| ). 



V. CONCLUSIONS 

In summary, following earlier studies!'! stimulated by Widom,! we have investigated 
the singular shapes of the various isothermal binodals, or two-phase coexistence curves, in 
the density plane near a critical endpoint. However, whereas the previous studies assumed 
classical or van der Waals expressions for the critical thermodynamics, our work is based 
on nonclassical phenomenological scaling postulates set out, in Sec. y, in a general form 
encompassing a spectrum of correction-to-scaling variables. Four types of critical endpoint s 
were distinguished and examined in detail: nonsymmetric^ labeled NA or NB depending on 
whether the lambda line Tx[g)^ which terminates at the endpoint [g^^ Tg), lies. A, below T = 
Tg (as in Fig. ||) or, B, runs above (as in Fig. H); and symmetric^ labeled, correspondingly, 
SA and SB: see Fig. ^. At the endpoints, the lambda-line binodals B^^ and — see 
Fig. |^(b) — were found to be singular with a leading "renormalized" exponent (1 — a)//? 
and subdominant singular correction exponents. The symmetric A binodals are displayed in 
( p.lOp [with explicit amplitude expressions recorded in ( |4.26| )-( ^y2^ )]; the nonsymmetric A 
binodals are presented in ( p.l7|) . 

Second, the noncritical or spectator-phase endpoint binodals B^^ and B'^~ — see 
Figs. 1(b) and |(b) — were found to be singular with a leading exponent {S + 1)/S (as 
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conjectured by Widomi): the symmetric binodals are given in (|3.36|) with the closely spaced 
sequence of correction exponents listed in (|3.38|) . The nonsymmetric binodals are similar 
but more complicated: see ( |3.44| )-( pI^9D . The endpoint binodals i3f and B2 which limit the 
critical phases — see Figs, ^(c) and |^(c) — have also been studied and found to have the 
same leading exponent, (1 — a)//?, as the lambda-line binodals: the symmetric forms are 
given in ( |3.52|) , the nonsymmetric expressions in (|3.56| ). 



In addition, above the endpoint temperature the binodals separating the spectator-phase 
from the near-critical phase — see and B'^t in Figs. |(c) and |(c) — were studied. 
They are analytic, but their slopes and curvatures develop singularities as T ^ Te+. The 
spectator-phase binodal is given in ( |3.58|) and (|3.59| ): its curvature diverges like (T — Tg)"^ 



when the critical endpoint is approached. The conjugate, near-critical-phase binodal is 
described by (|3.64| ) and ( p.67|) ; but in this case both the slope and the curvature vanish^ 
although in a singular fashion, on approaching the endpoint. 

Finally, the binodals that approach the three-phase region below the endpoint temper- 
ature have been considered. The spectator-phase binodals and — see Figs. |^(a) 
and|^(a) — are presented in ( |3.68|) and ( p. 69]) : as above the endpoint, their curvatures both 



diverge when T ^ Tg. 

Our analysis has utilized certain essential convexity or thermodynamic stability prop- 
erties at and near a critical endpoint and, for Ising-type criticality, also invoked a specific 
positivity of a scaling function expansion coefficient: see the discussion after Eqns ( |2.20|) 
and (|2.24| ). These features are taken up in Ref. 17. 
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FIGURES 

FIG. 1. The thermodynamic field space [g^T^h) showing a norzsymmetric (N) critical end- 
point, at the meet of a A line, marking the edge of a phase boundary surface p on which phases 
[3 and 7 can coexist, and a phase boundary surface a limiting the spectator phase a. The triple 
line r, on which phases a, (3 and 7 may coexist, extends above T = into the dot-dash line f 
which is the intersection of a with the extended phase boundary p (not shown). Note, as discussed 
below, that the A line shown here slopes downwards towards the a phase as T rises, so representing 
case A. 

FIG. 2. Isothermal sections of the NA endpoint phase diagram in the field space in Fig. 1 for 
(a) T < Te, (b) T = Te, and (c) T > T^. For T <Te the phase boundary a breaks into two pieces: 
(j+ separating phases a and and cr_ separating a and 7. The dotted curve represents the locus 
h{g^ T^h) = [see Sec. II] which coincides with the surface p in Fig. 1 and defines its extension p 
and, hence, the extended triple line f. 

FIG. 3. Isothermal density-density (or composition) diagrams for (a) T < Tg, (b) T = Tg, and 
(c) T > Te for an NA endpoint showing single-phase regions a, /3, 7, and /37, two-phase regions 
ruled by tie-lines connecting coexisting phases, and a three-phase triangle (dotted) in which coexist 
phases corresponding to the vertices ra, t/3, and T7. The various analytically distinct binodals 
are labeled Bj, • • where the superscript indicates the phase bounded by the binodal while 
the subscript serves (as needed) to specify the temperature, T ^ Tg. The same notations apply to 
a symmetric SA endpoint. At T = Tg the endpoint tie-line E^E^ defines the m or m = axis, 
shown dashed, where m and fh are fixed linear combinations of the densities pi and p2 (see Sec. 
Ill); the m and m axes on the plots in (a) and (c) have been omitted for the sake of clarity but 
are useful to understand the motion of the various features as T passes through Tg. 

FIG. 4. Thermodynamic field space illustrating a syrrmietiic critical endpoint, SB, for case B 
in which the A line slopes upwards away from the spectator phase a as T increases. Beyond these 
differences, the phases, phase boundaries, etc., correspond precisely to those in Fig. 1. 



47 



FIG. 5. Isothermal sections of the SB endpoint phase diagram in Fig. 4 for (a) T < Te, (b) 
T = Te^ and (c) T > T^. Compare with Fig. 2 and note that for T > the A point and its phase 
boundary p are disconnected from the boundary a. 

FIG. 6. Isothermal density-density diagrams for (a) T < Tg, (b) T = Tg, and (c) T > for 
the SB type of critical endpoint shown in Figs. 4 and 5. Compare the dispositions of the binodals 
with those shown in Fig. 3 and note the augmented labelling notation. 
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